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ABSTRACT 

Analyzing future weak lensing data sets from KIDS, DES, LSST, Euclid, WFIRST requires 
precise predictions for the weak lensing measures. In this paper we present a weak lens- 
ing prediction code based on the Coyote Universe emulator. The Coyote Universe emulator 
predicts the (non-linear) power spectrum of density fluctuations (Pg) to high accuracy for 
k G [0.002; 3.4] h/Mpc within the redshift interval z e [0; 1], outside this regime we extend 
Ps using a modified Halofit code. 

This pipeline is used to calculate various second-order cosmic shear statistics, e.g., shear 
power spectrum, shear-shear correlation function, ring statistics and COSEBIs (Complete Or- 
thogonal Set of EB-mode Integrals), and we examine how the upper limit in k (and z) to which 
Ps is known, impacts on these statistics. For example, we find that fcmax ~ 8 h/Mpc causes a 
bias in the shear power spectrum at{~ 4000 that is comparable to the statistical errors (in- 
trinsic shape-noise and cosmic variance) of a DES-like survey, whereas for LSST-like errors 
^max ~ 15 h/Mpc is needed to Umit the bias at £ ~ 4000. 

For the most recently developed second-order shear statistics, the COSEBIs, we find that 9 
modes can be calculated accurately knowing Pg to k^^ax = 10 h/Mpc. The COSEBIs allow for 
an EB-mode decomposition using a shear-shear correlation function measured over a finite 
range, thereby avoiding any EB-mode mixing due to finite survey size. We perform a de- 
tailed study in a 5 -dimensional parameter space in order to examine whether all cosmological 
information is captured by these 9 modes with the result that already 7-8 modes are sufficient. 
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1 INTRODUCTION 

Weak lensing by large-scale structure, so-called cosmic shear, was 
first detected in 2000 independently by several groups (Bacon et al. 
2000; Kaiser et al. 2000; van Waerbeke et al. 2000; Wittman et al. 
2000), and has recently progressed to an important tool in cos- 
mology. Latest results (e.g. van Waerbeke et al. 2005; Semboloni 
et al. 2006; Hoekstra et al. 2006; Schrabback et al. 2007; Hetter- 
scheidt et al. 2007; Massey et al. 2007b; Fu et al. 2008; Schrabback 
et al. 2010) already indicate its great ability to constrain cosmolog- 
ical parameters which will be enhanced by large upcoming ground 
based surveys like KIDS' , DES^, LSST^^, and satellite missions Eu- 
clid" andWFIRST^ 

In order to meet the requirements for analyzing high precision data 
sets from these future surveys cosmic shear needs to overcome sev- 
eral challenges. On the observational side accurate shape measure- 
ments (see Heymans et al. 2006; Massey et al. 2007a; Bridle et al. 
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2009; Bernstein 2010, for latest developments) and photometric 
redshift information (Bernstein & Huterer 2010; Hildebrandt et al. 
2010) probably pose the most challenging problems. On the astro- 
physical side intrinsic alignments can mimic a shear signal (Hi- 
rata & Seljak 2004; Mandelbaum et al. 2006; Joachimi et al. 2010; 
Mandelbaum et al. 2010) and need to be either removed (King & 
Schneider 2003; Joachimi & Schneider 2008, 2009), or modeled 
carefully in the subsequent analysis (Bridle & King 2007; Joachimi 
& Bridle 2009; Schneider & Bridle 2010; Kirk et al. 2010). 
An important step to quantify and check for these systematics is the 
decomposition into E- and B-modes, where, to leading order, grav- 
itational lensing only creates E-modes. In principle B-modes can 
arise from the limited validity of the Born approximation (Jain et al. 
2000; Hilbert et al. 2009) or redshift source clustering (Schneider 
et al. 2002). Predictions coming from numerical simulations differ 
on the impact of these effects (e.g. Heavens et al. 2000; Crittenden 
et al. 2001; Jing 2002), however, the observed B-mode amplitude 
is higher than one would expect from the foregoing explanations. 
Lensing bias (Schmidt et al. 2009) is another possible source of 
B-modes; again the amplitude compared to the E-mode is small 
(Krause & Hirata 2010). Most likely, a B-mode detection indicates 
remaining systematics in the observation and data analysis, in par- 
ticular an insufficient PSF-correction. 

There exist several methods to perform the EB-mode decomposi- 
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tion; as a general classification they can be separated into real- and 
Fourier-space methods, where the latter are inspired by the pseudo- 
CL formalism invented to analyze the CMB polarization power 
spectrum (Hivon et al. 2002; Brown et al. 2005). A finite survey 
area and masking eff^ects introduce a mixing of E- and B-modes 
in the CL's (also called leakage), which prohibits a clean separa- 
tion (Lewis et al. 2002). The effect cannot be removed completely, 
however for CMB polarization the arising leakage B-mode can be 
suppressed (Lewis 2003; Smith 2006; Kim & Naselsky 2010) to a 
level that enables for a possible detection of primordial B-modes 
(depending on the primordial tensor-to-scalar ratio). 
Hikage et al. (2010) extend the Pseudo-CL formalism to weak lens- 
ing, and test it on ray-tracing simulations, finding an B-mode leak- 
age at the percent level and below depending on the Fourier mode 
{. However, the strength of this effect depends on masking and sur- 
vey geometry and it needs to be examined on a case by case basis. 
Compared to CMB polarization additional difficulties arise in lens- 
ing when calculating the errors/covariances of the Pseudo CL's. For 
CMB polarization this covariance can be expressed in terms of the 
Pseudo CL's themselves assuming that the underlying field is Gaus- 
sian (Challinor & Chon 2005). This assvunption is justified for the 
CMB field (if primordial non-Gaussianity is small/zero), however 
it is not true for the late time shear field where non-linear structure 
growth leads to non-negligible non-Gaussian effects (White & Hu 
2000; Semboloni et al. 2007; Takada & Jain 2009; Sato et al. 2009). 
Here, higher order terms arise in the Pseudo CL shear covariance 
and the impact of EB-mode leakage on these terms is unknown. 
Several real space EB-decomposing methods sufi^er from EB-mode 
mixing as well, e.g., the aperture mass dispersion, the shear E- 
mode correlation function, and the shear dispersion (Kilbinger et al. 
2006). These statistics can be calculated in terms of the shear two- 
point correlation function (2PCF), however all three measures need 
information on scales outside the interval [<>min;'?max] of the mea- 
sured 2PCF. This information can of course be modeled using a 
theoretical 2PCF, however this biases the results towards the cos- 
mological model assumed in the 2PCF extension. The ring statis- 
tics (Schneider & Kilbinger 2007; Eifler et al. 2010; Fu & Kilbinger 
2010) and more recently the COSEBIs (Schneider et al. 2010, here- 
after SEKIO) provide a new method to perform an EB-mode de- 
composition using a 2PCF measured over a finite angular range, 
thereby avoiding any EB-mode leakage/mixing. 
Independent of whether cosmic shear data is analyzed in Fourier 
or real space both methods rely on accurate predictions for the cor- 
responding shear measure in order to constrain cosmological mod- 
els to the desired precision. Huterer & Takada (2005) find that one 
needs to know the power spectrum of density fluctuations {Ps) to 
<1% accuracy over a range of 1-10 h/Mpc in order to obtain cos- 
mic shear predictions that arc sufficiently accurate for LSST. The 
Halofit fitting formula for P,, described in Smith et al. (2003) is 
only accurate to ~10% (depending on cosmology and scale). The 
Coyote Universe emulator (Heitmann et al. 2009, 2010; Lawrence 
et al. 2010) has significantly improved on this issue; it models Ps 
with percent accuracy in a 5 dimensional parameter space over the 
range k e [0.002; 3.4] h/Mpc within z € [0; 1]. 
This paper has two goals: First, we built a predictions pipeline for 
second-order shear measures based on the Coyote Universe emula- 
tor and a modified Halofit. For improved weak lensing predictions 
from future numerical simulations we examine to which scale k and 
redshift z the density power spectrum must be known to model var- 
ious second-order shear statistics to the desired accuracy for DES 
and LSST. The second goal is to examine the COSEBIs in a 5- 
dimensional parameter space. As outlined in SEKIO the COSEBIS 



do not only offer a check for B-mode systematics, moreover the au- 
thors argue that the COSEBIs' E-mode is the quantity that should 
enter a cosmic shear likelihood analysis. The 2PCF itself in par- 
ticular should not be used for this purpose. The reason for this is 
that even when finding no B-modes in any EB-mode decomposi- 
tion tests, the 2PCF can still be affected by B-modes that 1) mimic 
a constant or linear shear field 2) are localized in the power spec- 
trum but present on scales {€) larger than the survey area. 
The (logarithmic) COSEBIs are designed to contain all EB- 
decomposable cosmological information, and compress it into few 
data points (modes). SEKIO examined this in the two-dimensional 
parameter space fl^ vs. cr^; we extend their analysis to 5 parame- 
ters which correspond to the cosmological parameter space covered 
by the Coyote Universe. 

The paper starts with a short description of our weak lensing pre- 
diction pipeline; here and in Sect. 3, we also determine the require- 
ments for future numerical simulations needed for sufficiently ac- 
curate weak lensing predictions for DES and LSST. In Sect. 4 we 
examine the performance of the COSEBIs in a 5-dimensional pa- 
rameter space, in particular we address the questions 1) how many 
modes of the COSEBIs need to be included before the cosmologi- 
cal information saturates and 2) how does this saturation limit com- 
pare to the case of a pure E-mode 2PCF (which we can of course 
simulate but is never guaranteed in a real data set). We discuss our 
findings and conclude in Sect. 5. 



2 THE SECOND-ORDER SHEAR STATISTICS 
PREDICTION CODE 

The code is based on the Coyote Universe power spectrum emu- 
lator (Lawrence et al. 2010), which in its newest version emulates 
Ps over the range k e [0.002; 3.4] h/Mpc within z e [0; 1]. For 
the higher k and z range we calculate the linear power spectrum 
from an initial Harrison-Zeldovich power spectrum (Ps(k) oc kf^) 
using the transfer function of Efstathiou et al. (1992). We account 
for the non-linear structure growth using the formula of Smith et al. 
(2003) and then scale the result such that the transition to the Coy- 
ote Universe power spectrum is smooth. More precisely, we derive 
the scaUng factor by taking the ratio of (Coy ote)//'^ (Halofit) for 
the largest k of the Coyote emulator. For all k > femax (Coyote) we 
multiply the amplitude of Halofit with this factor. 
The main difficulty when co-implementing the emulator and 
Halofit is that the former automatically calculates Hq such that it 
fits the WMAP-5year constraints on the angular scale of the acous- 
tic peak (Komatsu et al. 2009), which implies that our code cannot 
vary f2„, fit independently of Hq. 



2.1 The shear power spectrum 

The shear power spectrum P^ can be expressed as a an integral over 
the density power spectrum Ps through Limber's equation 

4c^ Jo a\w) \Mw) I 

where { is the 2D wave vector perpendicular to the Une of sight, w 
denotes the comoving coordinate, Wh the comoving coordinate of 
the horizon, a(w) is the scale factor, fKiw) the comoving angular 
diameter distance. 

The weight factor g is defined as an integral over the redshift dis- 
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Figure 1. The shear power spectrum when calculated with Halofit or from 
the new Coyote Universe/modified Halofit code. 



tribution of source galaxies n{w(z)) 

g(w) = I aw n{w )- 



w) 



In the calculation of we choose a redshift distribution of source 
galaxies similar to that for the CFHTLS described in Benjamin 
et al. (2007) 



n{z) 



— exp 



zo 



(3) 



withtt = 1.197,j8 = 1.193, zo = 0.555, and a cut-off at z^^x = 4. 
Figure 1 compares the shear power spectrum calculated from 
our modified code compared to using Halofit only. The as- 
sumed cosmology corresponds to the best fit cosmology 
of the WMAP 7-year analysis (Komatsu et al. 2010), i.e. 

= 0.272, = 0.728, erg = 0.807, h = 0.703, fib = 0.045, 
«s = 0.961, Wo = -1, and w„ = 0. These parameters also define 
our fiducial cosmology used for the fiducial data vector in the 
likelihood analysis (Sect. 4). A direct comparison to power spectra 
from raytracing simulations is unfortunately not possible due to 
the fact that fi„,, fit. cannot be chosen independently, hence there 
was no raytracing simulation with the "correct" cosmology at hand. 



2.2 Required accuracy in Ps for future weak lensing surveys 

For a future extension of the Coyote Universe project it is valuable 
to determine exactly over which scales k and redshifts z the density 
power spectrum must be known to the 1% accuracy requirement 
defined in Huterer & Takada (2005). Obviously the redshift range 
will depend on the redshift distribution of the considered survey; 
we outlined our choice in Sect. 2.1 and argue that this is likely to 
be close to what is expected for DES. In Fig. 2 we plot the ratio 
Pe{z)IPe(z = 4) as a function of z, where Pe(z) is calculated as in 
Eq. 1 but with w(z) as the upper limit. The shaded regions indicate 
statistical errors (intrinsic shape noise and cosmic variance) 



1 



Pe + 



2n„,, 



gal 



(4) 



for aDES-like survey (light shaded outer region, /,ijy=0.125, Wgai = 
I2/arcmin^), and a LSST/Euclid-like survey (dark shaded inner re- 



Figure 2. The fraction of Pe(z)/Pe{z = 4) for various different C where 
Pe(z) is calculated from Limber's equation but with w(z) being the upper 
limit of the integral. We see that the shear power spectrum can be calculated 
accurately when knowing Ps only up to z = 1.5. This will of course depend 
on the assume redshift distribution. The error bars are calculated for a DES- 
like survey (light shaded), and for a LSST-like survey (dark shaded) atl = 
1000. 



Table 1. The f in Pe where the eiTor in the shear power spectrum A = 
l^'E(^max)/^'E(^max = 100) - 1| exceeds a certain percent level (see also Fig 
3). The values of k are given in h/Mpc. 



A 


^max — 15 


^max — 10 






1 % 


e = 3593 


{ = 2096 


f = 1557 


( = 311 


5 % 


e = 6604 


^ = 4152 


f = 3152 


e= 823 


10% 


e = 8589 


e = 5483 


£ = 4213 


^= 1162 



gion, /sky=0.5, Wgai = 40/arcmin^). We choose cTc = 0.3 in both 
cases. 

Here, the error bars are calculated at f ~ 1000 for the reason that 
they become minimal around this scale (see Fig. 3). We find that 
for the considered redshift distribution, it is sufficient to know P^ 
to z ~ 1.5 even for LSST like error bars. We note however that due 
to the increased depth of an LSST like survey its redshift distribu- 
tion might be significantly different from what we assumed here. 
Next we impose a cut-off scale k^.^^ in P^, i.e. we set P^ = 
for k > /:n,ax, and study the impact of different choices of A;max 
on the shear power spectrum and the shear 2PCF. The top panel 
in Fig. 3 shows the relative error in the shear power spectrum 
APe = |/'E(^max)//'E(fcmax = 100)| for various k^.^^. We find that 
for DES-like error bars A;,^^^ = 8h/Mpc is sufficient if one is only 
interested in P^it) below C = 3960, fc^ax = lOh/Mpc extends this 
range to{ = 5995, and if Ps can be modeled up to A'„mx = 15h/Mpc 
even £ = 15000 is still within the DES error bars. This changes con- 
siderably for LSST, where knowing Ps only up to ky„,,^ = 8h/Mpc 
gives uncertainties in Pe(^) above € = 1750 that are equal/larger 
than the LSST errors. For ky„.^^ = lOh/Mpc this threshold is shifted 
to f > 2308, and for /t.^^^ = 15h/Mpc to { > 4127. 
We perform a similar analysis for f+ (Fig. 3, middle panel) and f_ 
(Fig. 3, lower panel) which are calculated from P^ via 

f±(i?)= r ^h/4m[PE(n + PBm ■ (5) 



The corresponding error bars for f+ are derived from diagonal ele- 



ments of the shear covariance A^ = Je, which can be calculated 
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Figure 3. Upper: The relative error in the shear power spectrum 
^'E(^max)/^'E(^max = 100) for various cut-off scales ^niax in the density 
povi'er spectrum. The light shaded region indicates the statistical errors as- 
sumed for a DES-like survey, the dark shaded region conesponds to LSST- 
like en'or bars. Middle: The 2PCF f + and Bottom:^- (right) for various /:,^iix 
in Ps- Here we only show DES-like errors. 



in terms of the power spectrum (Joachimi et al. 2008) 

Cl,0?„)?,) = r d^no/4(^J?,)Jo/4(^'?;)f/'EW + ^ 

(6) 

The above expression only holds under the assumption that the four 
point function of the shear can be expressed in terms of power 
spectra which again only holds for Gaussian shear fields; including 
non-Gaussianity increases these error bars (Semboloni et al. 2007; 



Takada & Jain 2009). More importantly, the 2PCF covariance has 
significant off-diagonal terms, hence its errors are strongly corre- 
lated which implies that the ratio of error bar to signal, as shown in 
Fig. 3, does not fully resemble the actual signal-to-noise ratio of the 
2PCF. For this one has to perform a full likelihood analysis taking 
the total covariance into account (see Sect. 4). 
Nevertheless our results in Fig. 3 show that f + is hardly affected by 
a cut-off even when choosing fcniax = 3h/Mpc. In contrast the im- 
pact on especially on small scales, is much more severe. Even 
^max = 15 h/Mpc is hardly within the DES error bars on scales < 1'. 
The insensitivity of f+ to the choice of kmax can be explained when 
looking Eq. 5 and the filter function Jo(^i?). We see that the main 
contribution to the integral in Eq. 5 comes from small £ (where Pe is 
not affected by a cut-off in Ps), and that this tendency increases for 
larger arguments of Jq, i.e. larger i?. In contrast, for f_ we see that 
]4({&) strongly suppresses these small i and collects power mainly 
from larger scales, which is the reason why it is more biased when 
using low 



3 COSEBIS AND RING STATISTICS 

We extend the analysis to the most recently developed EB- 
decomposition methods, namely the ring statistics and the COSE- 
BIs. Both can be expressed as integrals over the 2PCF, i.e. 



-'•'mi,, 

for the ring statistics, and 



r'™'" dj? r 

1„, 



.T'°^{&)U-ff)\ 



(7) 



(8) 



in case of the COSEBIs. The filter functions Z± are derived in 
Schneider & Kilbinger (2007), we use the modified version de- 
scribed in Eifler et al. (2010). The derivation of the COSEBIs' fil- 
ter functions T+„ is outlined in detail in SEKIO. The superscript '°= 
indicates that the roots of the corresponding filter function are dis- 
tributed logarithmically in ??. We note that a similar relation exists 
for linear T-functions, which can be expressed very conveniently in 
terms of Legendre polynomials. However T'l'Jf has the advantage to 
compress the information contained in the 2PCF into significantly 
fewer data points hence we consider only r+°f in this paper. 
Another way to calculate the COSEBIs EB-mode employs the 
shear E-mode power spectrum, or a B-mode power spectrum, re- 
spectively 



At I 



E„IB„= --p^^^(^[)w':He) . 



(9) 



where the filter functions w]?^ can be calculated from the ri°f as 



(10) 

(11) 



A similar relation exists for the ring statistics but using Z+ instead 

Similarly to Eqs. (8, 9) the covariance of the COSEBIs' (and ring 
statistics') E-mode can be expressed through the covariance of the 



© 2010 RAS, MNRAS 000, 1-9 



Weak lensing statistics from the Coyote Universe 5 



and the covariance of Pe 

1 I cP- \^ 

47r^/sky Jo \ 2«gal/ 

= - I d!? 1? I di>' 0' 

X 2 T'^mT'Z\&')C,M,-»') . (13) 

In Fig. 4 we perform a similar analysis as in Fig. 3 but for the ring 
statistics {left panel) and for the COSEBIs {right panel). The error 
bars for the ring statistics are calculated as the square root of the di- 
agonal elements of the ring statistics covariance. In contrast to the 
2PCF this gives a fair estimate of the signal-to noise ratio, as the 
ring statistics covariance is almost diagonal (Eifier et al. 2010). We 
find only a small bias depending on kmx\ only for kmsn = 3 h/Mpc 
there is a clear deviation, all other estimates of the ring statistics 
signal are well within the DES error bars. This behavior can be ex- 
plained again when looking at the ring statistics as a filtered version 
of the shear power spectrum, and the corresponding filter function 
W{f) (Schneider & Kilbinger 2007, Fig. 3). The filter function sup- 
presses the large { information, where the shear power spectrum is 
most affected by a low fcmax- 

We perform a similar analysis for the COSEBIs, however refrain 
from plotting the error bars as they are heavily correlated, at least 
for the logarithmic case. We postpone an error analysis for the 
COSEBIs to the next section, where we carry out a complete like- 
lihood analysis in a 5-dimensional cosmological parameter space, 
that takes the full COSEBIs covariance into account. Nevertheless, 
one can see that the COSEBIs can be calculated accurately up to the 
9th mode for fe^ax = 10 h/Mpc. Choosing fcmax = 3 h/Mpc results in 
a notable difference from the 7th mode on. It is the major purpose 
of the next section to check how many modes are needed to capture 
the cosmological information of second-order shear measures in a 
5-dimensional parameter space. 



4 LIKELIHOOD ANALYSIS 

In SEKIO we analyze the performance of the COSEBIs relative 
to the theoretical 2PCF for the Slm-crg parameter space. Here, we 
extend this analysis in several ways: First, we consider a parame- 
ter space consisting of 5 parameters (w„, wj,, erg, n^, wo) to quantify 
how many modes Sjf* are needed to capture the bulk of cosmo- 
logical information in this increased parameter space. Note that we 
use a)„, = Si^h^ {<JJh = Sibh^) instead of fi,^ {^b), for the reason 
explained in Sect. 2. Second, we compare the information con- 
tent of the COSEBIs when calculated from a 2PCF that is mea- 
sured on only small scales (i? 6 [1'; 100]) and only large scales 
{& e [120'; 400']). This analysis is motivated by the findings in 
SEKIO (see Fig. 11), namely that for fi^im = 1' the information 
content of the COSEBIs does not increase significantly when in- 
creasing i?max beyond 200' . 

The data vectors of the COSEBIs are calculated from the 2PCF 
via Eq. (8), and we cross-check the results using Eq. (9). Similarly 
we calculate the covariance using both Eqs. (13) and (12). For co- 
variances this consistency check is more important as numerical 
difficulties arise more likely in the covariance calculation (see Ap- 
pendix of SEKIO for details). The survey parameters we assume in 
the calculation of all covariances in this paper correspond to those 
of the Dark Energy Survey (see Sect. 2). 




Figure 5. Tfie dependence of q and / on the survey size for the Hm-crg pa- 
rameter space only. This illustrates the difference of a fisher matrix based 
figure-of-merit and a figure-of-merit that is based on a full likelihood anal- 
ysis. 

4.1 Likelihood formalism 

In this paper we quantify the information content of a measure us- 
ing the figure-of-merit 

q=^/\Q\ with Q,j^ jSnp{n,j\d){7T,-K]){nj-ji]). (14) 

In our case 2 is a 5 X 5 matrix, which can be interpreted as the ex- 
pectation value of the parameter covariance matrix. The parameter 

vector reads ti = {co,,,, cot, n,, erg, wq), and n^. and n^j denote the fidu- 
cial parameter values. Note that the posterior likelihood p{nij\d) is 
calculated by marginalizing over the likelihood of the three param- 
eters nuMj. 

The figure-of-merit q corresponds to the more common Fisher ma- 
trix based figure of merit / = 1/| VF| (see Tegmark et al. 1997, 
SEKIO for the exact definition) if the likelihood in parameter space 
is Gaussian. The Fisher matrix F can be interpreted as the expec- 
tation value of the inverse parameter covariance evaluated at the 
maximum likelihood estimate parameter set, which in our ansatz 
corresponds to the fiducial parameters. Mathematically we can ex- 
press this equivalence as 

f=-|^=^Jic:\=M = q■ (15) 
ViFf 

However, the assumption of a Gaussian likelihood in parameter 
space is only justified close to the maximum of the likelihood func- 
tion, i.e. close to the fiducial parameters. If the likelihood func- 
tion does not fall off quickly enough, / is a rather bad approxi- 
mation and can give significantly different values compared to q 
(see SEKIO). Obviously the Gaussian likelihood assumption signif- 
icantly becomes more valid the larger the survey volume and we ex- 
amine this behavior quantitatively in Fig. 5 for the two-dimensional 
parameter space n„,-cr8. We find that q and / deviate significantly 
for small survey volumes, whereas this difference drops below 1% 
for surveys > 9300 deg^. Here, we consider a survey size corre- 
sponding to the Dark Energy Survey, i.e. 5000deg^ where the devi- 
ation is still >7 %; therefore we refrain from using / in our analysis 
and consider q only. 

In the calculation of the posterior likelihood we assume that the 
data points in d (which in our case are either COSEBIs or 2PCFs) 
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Figure 4. The ring statistics {lejf) witli DES-errors and the COSEBIs {right) for various /cniax in Pg- Note that we refrain from plotting errors for the COSEBIs 
as they are highly correlated (see text). 



are following a Gaussian error distribution 



1 



(2;r)"/2 VlC| 



exp 



(16) 



and that the covariance C is constant in parameter space. These 
assumptions are not fully justified (see Hartlap et al. (2009); Seo 
et al. (2010) for information on the first and Eifler et al. (2009) 
for information on the second assumption) and they need to be 
examined carefully when inferring cosmological constraints from 
data. For our purpose these assumptions should not influence the 
results qualitatively as we are comparing the information content 
of 2PCF and COSEBIs relative to each other. 



4.2 Results 

In Sect. 4.1 we describe our figure of merit as the determinant of 
the parameter covariance matrix Q. In addition, to the full parame- 
ter covariance we calculate the determinants of several submatrices 
of Q; these q can be illustrated as two-dimensional likelihood con- 
tours where we marginalize over the other parameters. 
In Fig. 6, we show the q of the COSEBIs compared to the pure 
E-mode 2PCF for the 5-dimensional parameter space {top left) and 
when considering the subsets. w„ vs. wo {top right), (o^ vs. erg {bot- 
tom left), M, vs. Wo {bottom right). We performed the same analysis 
for the other 6 parameter subsets, but only show this selection due 
to the fact that the results look very similar. Note that scale of q 
is logarithmic in the top left panel whereas it is linear in all other 
cases. In addition, we scale all q such that ? = 1 for the [r;400']- 
2PCF(dashed line). 

In each panel we consider three intervals in angular scale [1'; 400'], 
[1'; 100'], and [120'; 400'], and calculate the value of q for the E- 
mode 2PCF as the upper limit on cosmological information from 
second-order shear measures. Note that the line representing the 
[1'; 100']-2PCF overlaps with the line of the [l';400']-2PCF We 
then calculate the q of the COSEBIs as a function of modes that 
are included in the likelihood analysis, and thereby quantify their 
cumulative cosmological information as a function of modes. 
Considering the [1'; 400'] -interval we find that even in a 5- 
dimensional parameter space the number of COSEBI-modes that 
is required to capture the cosmological information is quite small, 
~7-8 modes {top left). This guarantees a manageable covariance 



matrix and a stable calculation of its inverse. In addition, we see 
that knowing the density power spectrum to = 10 h/Mpc is suf- 
ficient to avoid uncertainties in the data analysis coming from the- 
oretical predictions. The difference of the 8-mode COSEBIs to the 
shear 2PCF is very small, hence the loss in information is a negli- 
gible inconvenience compared to the uncertainties that are avoided 
when using the COSEBIs instead of the 2PCF. 
Slightly surprising is the difference in information when com- 
paring the [1'; 100']- and the [120';400']-interval. The [1'; 100']- 
COSEBIs' information content saturates at ~5 modes for 5 cos- 
mological parameters, and the difference to the [1'; 400']-case is 
rather small. Comparing the q of the [1'; 100']- and [1'; 400']-2PCF 
we can hardly see a difference, which indicates that the bulk of 
cosmological information is indeed contained on small scales. The 
[120';400']-COSEBIs saturate already at 2-3 modes, and the dif- 
ference to the other saturation limits is large. Also the difference 
between COSEBIs and 2PCF is much more significant in this case 
compared to the other two. 

To explain this behavior one needs to examine the variation of the 
signal-to-noise ratio of the 2PCF with respect to cosmology. Know- 
ing that the bulk of the cosmological information is contained in 
and assuming that the square root of the diagonal of the covariance 
serves as a rough estimate for the noise-level, we see from Fig. 3 
that the signal-to-noise is relatively low on scales >100' compared 
to the scales <100'. Plotting the variation of the 2PCF signal in 
the 5-dimensional cosmological parameter space we find the rel- 
ative deviation to the signal of the fiducial cosmology is roughly 
constant over all scales. The lower noise level on small scales then 
impUes that these scales contribute more information on cosmol- 
ogy- 
There are only slight variations to this behavior when looking at the 
other panels of Fig. 6. As expected the information content of the 
COSEBIs saturates earlier compared to the 5-dimensional parame- 
ter case, however the fact that almost all the information is present 
on scales <100' remains. 



5 CONCLUSIONS 

Accurate predictions of cosmic shear measures play an important 
role in the analysis of future weak lensing data sets. In the first part 
of this paper we present our new weak lensing predictions pipeUne 
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Figure 6. The information content of the 2PCF and the COSEBIs for the 5-dimensional parameter space (top left) and for three two-dimensional parameter 
spaces quantified by the measure q. Note that smaller ^corresponds to larger information content (smaller hkeUhood contoitts). The lines correspond to the 
2PCFs' q for the 3 considered !?-intervals and pose the upper limit of the second order information content. The values of q for the COSEBIs are plotted as a 
function of number of modes n, that are included in the likelihood analysis. 



that is based on the Coyote Universe emulator and extends the emu- 
lator using the Halofit code. A shear power spectnmi derived from 
this new pipeline differs from a shear power spectrum calculated 
from Halofit only by 6-11 %, depending on the Fourier mode (. 
We consider this as a first step in developing a weak lensing pre- 
dictions pipeline that meets the requirements for DES and later for 
LSST, Euclid, and WFIRST. For this we examine to which k in 
h/Mpc and z, the density power spectrum must be modeled ac- 
curately by numerical simulations in order to meet those require- 
ments. We find that A:max = 8 h/Mpc causes a bias in the shear 
power spectrum at £ = 4000 that is within the statistical errors 
(intrinsic shape-noise and cosmic variance) of a DES-like survey, 
whereas for LSST already femax = 15 h/Mpc is needed. 
A future pipeline for weak lensing predictions that models the den- 
sity power spectrum to such large k, the shear power spectrum to 
such high {, respectively, will have to take baryons into account 
(Jing et al. 2006; Rudd et al. 2008). For example, Jing et al. (2006) 
find that the shear power spectrum on scales of ^ 6 [1000, 10000] 
are affected by 1-10%, depending on the level of complexity in the 
treatment of baryons. 

In addition to the treatment of baryons additional corrections need 

to taken into account. For example, shape distortions probe the re- 
duced shear^ = 7/(l-/<:) instead of the shear (y) itself (e.g., Shapiro 
2009), and when calculating the reduced shear from the measured 
mean ellipticity higher order terms in g need to be taken into ac- 



count. Corrections to the Bom approximation are necessary to ac- 
count for multiple deflections of light rays (e.g., Hilbert et al. 2009), 
and corrections to biases that occur because de(magmflcation) of 
galaxies due to lensing correlates with selection criteria of the con- 
sidered galaxy sample (e.g., Schmidt et al. 2009). These corrections 
are calculated in Krause & Hirata (2010); the authors find that re- 
duced shear and magnification bias are important already for DES- 
like surveys, whereas the other effects will only become important 
for weak lensing data from LSST. 

In the second part of this paper we extend earlier studies of the 
COSEBIs, the most recently developed EB-mode decomposing 
second-order cosmic shear measure. In particular, we are interested 
in their performance in a high-dimensional cosmological parameter 
space. Compared to other second-order shear statistics the COSE- 
BIs can be calculated from a shear 2PCF measured on a finite in- 
terval [i?mm; i^max]. The 2PCF again is independent of any masking 
effects or survey geometry and thereby avoids several difficulties 
present in other second-order shear measurement methods, e.g. the 
shear power spectrum. 

The COSEBIs can be imagined as a filtered version of the 2PCF 
(or the shear power spectrum); instead of having angular scale ■& or 

Fourier modes £ as an argument, they are a function of the order of 
the polynomial which is used as their filter function. Furthermore, 
they are designed to condense the second-order cosmic shear infor- 
mation into a (small) number of discrete modes. In SEKIO two dif- 
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ferent filter functions are examined: Filter functions that sample the 
2PCF linearly can be expressed very conveniently in terms of Leg- 
endre polynomials, however the number of modes that is needed to 
capture all the cosmological information is rather large in this case. 
The logarithmic filter functions are much more efficient; SEKIO 
find that ~ 5 modes are sufficient to capture the cosmological infor- 
mation in the two-dimensional parameter space Slra vs. erg. 
Here, we extend their analysis and examine the performance of 
the COSEBIs in a 5-dimensional parameter space. We find that 
the bulk of the cosmological information is contained on angular 
scale <100', whereas scales above this threshold only contribute 
little to constraints on cosmology. Furthermore, we find that ~8 
modes of the COSEBIs are sufiicient to capture all the EB-mode 
decomposable cosmological information. This number is still rel- 
atively small, in particular as we show that the COSEBIs can be 
calculated accurately up to the 9th mode when knowing Ps "only" 
up to krazx = 10 h/Mpc. From these results we conclude that the 
COSEBIs are fairly robust against theoretical uncertainties in mod- 
eling the density power spectrum, and represent a good choice as a 
second-order weak lensing statistics. 
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